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Eigenstates of the FQHE hamiltonian problem after to be projected on the LLL are determined 
for filling factors v — 1/q, with q an odd number. The solutions are found for an infinite class of 
finite samples in which the Coulomb potential is periodically extended. Therefore, a thermodynamic 
limit solution is also identified. The results suggest the presence of integrability properties in FQHE 
systems. The many particle states are simple Slater determinants constructed with special single 
particle states. These orbitals are defined as powers of order q of composite fermion like wavefunc- 
tions associated to a reduced magnetic field B/q. At the same time, those composite fermion states 
were obtained by factorizing and canceling fixed position (quasi-momentum independent) zeros in 
previously derived exact Hartree-Fock orbitals. A formula for the energy per particle of the FQHE 
states is given for finite samples as well as for the thermodynamic limit state. As a side result, the 
same composite fermions like orbitals are employed to construct variational wavefunctions of the 
system, showing zeros of order q as two electrons approach each other, as Laughlin states do. The 
long range spatial correlation associated to the starting HF solutions may further reduce the energy 
of these states. 

PACS numbers: 73.43.Cd, 73.43.-f 



I. INTRODUCTION 



The Fractional Quantum Hall effect has been a central theme of investigation in condensed matter physics for 
decadesi2^£&£. Up to now its understanding rests on Laughlin's variational approach for Landau level filling 
fractions of the form 1/q, for q odd, and its generalization to general fractional fillings developed by Jainii. The wave 
functions proposed by these authors are not eigenstates of the system hamiltonian although their energy is quite close 
to the true ground state, as numerical calculations have showr£. An alternative approach to FQHE originated in the 
special Hartree-Fock (HF) solutions proposed by F. Clare 8 ' 9 many years ago, which although lacking correlations, 
capture the essential odd-denominator rule observed in the fractional quantum Hall regime. Later, further work 
(Refs.— i 12 i 13 ' 14 i 15 ' 16 ) has shown promising properties of such solutions, raising the hope that by adding electron 
correlations one could approach the true ground state of the system, possibly leading to the long-standing idea of a 
weak dynamical breaking of the translation symmetry ^ 2 6 1 2 7 . 

In this work, in first place, we construct exact eigenstates of the restricted to the lowest Landau level (LLL) FQHE 
problem when the Coulomb potential is periodically extended. These wavefunctions corresponds to filling factors of 
the form v = 1/q, with q an odd integer. They satisfy periodic boundary conditions, and the Coulomb potential is 
defined to have the same periodicity. The possibility for the determination of these exact wavefunctions is offered by 
the use of a special basis of one particle states which were constructed by factorizing zeros in certain solutions of the 
HF problem. Such kind of functions were initially discussed in Ref . n ' 12 i 13 ' 14 i 15 ' 16 . The starting HF orbitals showed 
a structure, in which, from the total number of zeros (being equal to the number of quantum fluxes piercing the 
periodicity region) only a fraction q had positions depending of the quasimomenta quantum numbers of the orbitals. 
The rest of (q — l)/q zeros were situated at the origin of coordinates. Those fixed zeros where able to be factorized 
and canceled, by then leaving as a result a set of composite fermion like orbitals, associated to a reduced magnetic 
field B/q. This property directly led to the idea of the possible physical relevance of special single particle orbitals 
constructed as powers of order q of such states (to be called below: q product states). It can be noted that those 
products regain the number of zeros required by the actual magnetic field value. Further motivation of the line of 
search under consideration came after checking that the set of orbitals constructed for all the quasimomenta values, 
became closed under the action of the magnetic translations defining the periodically extended Coulomb potential. 
The extended potential is defined to have periods coinciding with the same assumed here von Karman periodicity 
lengths of the system: L\/q and L%/q. The mentioned closure property, then allows to directly show that the 
determinant functions formed with the q product states are in fact exact solutions of the Schrodinger equation. A 
simple formula for the exact eigen-energy of the system also follows. Afterwards, since the solution was determined 
for arbitrary but finite size of the systems, the thermodynamic limit of problem can also be determined. 
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The decision about whether these exact eigenstates of the FQHE systems are excited or the ground state ones, 
will be further considered in a work now in progress (paper II). That is, we do not argue that the presented exact 
solutions represent the ground state of the FQHE system, which the variational Jain-Laughlin approaches and the 
exact diagonalization results very well approximate. Occasionally, the determined states could result to be excited 
ones, showing even higher energy per particle than the Laughlin wavefunctions. However, even in that case, the 
wavefunctions could describe alternative phases of the 2DEG which having physical interest. This can be the situation, 
because the experiments in FQHE systems are performed over a wide range of conditions. By example, particularly 
important kind of measurements are related with finite temperature experience o 17 ' 18 i 19 ' 20 . In them, the existence of 
phase transitions between states of the 2DEG when the temperature changes has been detected. Such transitions are 
related with the relaxation of the FQHE systems to states which can not be described by just few excitation over 
the ground state, but should correspond two radically different ground states, showing higher energy per particle 
values. This circumstance indicates the physical interest of the investigated exact eigenstates after taken in the 
thermodynamic limit. On another hand, if the exact solutions could occasionally correspond to the ground states, 
their physical interest will be clear. 

It should be noticed that this initial part of the work contains a large amount of formal and technical developments 
which are however, required for arguing the exact nature of the wavefunctions. Therefore, in order to avoid a 
disproportionate larger paper, we plan to publish the results as a series of two articles: papers I and II. The present 
one (to be referred as paper I) will contain the formal derivation of the exact eigenstates. In the second work of the 
series (paper II) we plan to present the evaluations of the energy per particle and pair correlation functions of the here 
discussed eigenstates. Such results should define whether those wavefunctions correspond to phases that only can be 
realized at high temperatures, or on another hand, can also describe the exact ground states of the FQHE problem. 

It can be remarked that the connection of the exact states discussed here with the solutions of the exactly solvable 
model defined in a thin torus derived in Refs i 21 ' 22 i 23 ' 24 i 25 is also planned to be examined in the next paper II. 

In connection with the structure of the zeros of the obtained solutions, it can be underlined that the almost exact 
states determined in Refs . 32 i 33 ' 34 for Quantum Dots in FQHE configurations, indicate that the zeros of the exact 
wavefunctions are not of order q as it occurs in the Laughlin and Jain variational states. The determinant nature of 
the functions studied here, seems to be compatible with this property. However, this point is expected to be more 
clarified also in the planned paper II. 

Another objective of this work is treated in an Appendix. There, the same special one particle functions are also 
employed to construct ansatz for low energy eigenfunctions. They are defined as powers of order q of determinant 
like states associated to samples having a magnetic field reduced to a fraction q of the real external one. The powers 
q entering in the construction, ensure that the proposed states show the required number of zeros associated to the 
external field in the periodicity region. This form of definition directly implies that the zeros of the wavefunction at 
nearly coinciding coordinates Zi and Zj are of the form [zt — Zj) q . Thus, the advanced states show a similar short 
range behavior as the Laughlin states which thanks to such a property have the status of being the best variational 
state proposals up to now. However, since the starting point for the definition of the here advanced ansatz, are the HF 
one particle states, the expectation emerges, about that the long range correlations incorporated in the HF solutions, 
will remains to be effective in the new wavefunctions, by occasionally helping to lower the energies slightly more. 

The paper will proceed as follows. In Section II the wavefunction corresponding to an exact eigenvalue of the 
periodically extended Coulomb Hamiltonian reduced to the LLL, will be presented. Also, a formula for the energy 
per particle of the FQHE state is derived and the thermodynamic limit of the state is analyzed. Section III exposes the 
conclusions and a summary of the results. In Appendices A and B the definitions and auxiliary technical derivations 
are given for the sake of completeness and bookkeeping purposes. The construction of the ansatz wavefunctions for 
the ground state of FQHE systems at filling factors 1/q is discussed in Appendix C. 

II. EXACT SOLUTIONS OF THE FQHE PROBLEM 

Let us consider a 2-dimensional electron gas (2DEG)in the presence of a perpendicular magnetic field B such that 
a fraction v = 1/q of the lowest Landau level (LLL) is filled, where q is an odd integer. The electronic systems will 
be periodically extended in what follows to satisfy magnetic periodic von Karman conditions in the lattice generated 
by the vectors 4 - 

Li TV ai L2 N &2 

q q ' q q 

Here N is an even multiple of q, such that precisely one flux quanta pierces the cell of sides ~L\/{qN) and L12/N 
The unit cell vectors a! and a 2 define the smaller spatial magnetic translations under which the HF orbitals motivating 
this work, were eigenfuntion G n i 12 ' 13 i 14 ' 15 i 16 (see Appendix A for further definitions and conventions). 



3 



As it was mentioned in the Introduction, the Coulomb potential between particles will be periodically extended to 
be invariant under translations in the same lattice defining the von Karman periodicity of the electron wavefunctions 

TVai Na 2 

R-q = Rq.l h Rq,2 

q q 

— Rq.l ~ + Rq,2 — ) (1) 

q q 

in which R q .\ and R q ,2 are arbitrary integers. Then, the extended Coulomb potential has the expression 

^(x-x 0=£ .— ^TWJ . (2) 



|x-x' + R„ 



which is periodic in the lattice R q . After projection on the LLL, the Coulomb Hamiltonian can be written in the 
Landau gauge, in the form (See Appendix A and Ref. 



rEE^x-p(-^ T(i, z *)T(-i, z *), ( 3) 



, 2 7 re 2 r 2 , 

i i<j 

A c = ±^, (4) 
q 2 



where the magnetic translation operators in an arbitrary vector v are given by (See Appendix A) 



T(v,z*) = exp(-^ i (^+ ? z*))exp(^A), ( :„ 
v t 2 oz* 



and the vectors 1 arc defined in terms of the normal vector n and the discrete reciprocal lattice momenta (associated 
to the von Karman lattice) k as 

1 = r 2 n x k, (6) 
k = Tv s i + Tv s 2' ( 7 ) 

1 1 

q 1} q 2 = -oo, ..,-2,-1,0, 1,2, ...,oo, (8) 

in which the vectors Si and S2, are defined in Appendix A. Technical definitions and a derivation of the expression of 
the above periodically extended Hamiltonian can be also found in Appendix A. 

Let us now consider the following Slater wavefunction for a system of N e — (^) 2 electrons, 

^>(z* 1 ,z^...z Ne )=Det[Q li (z*)}, (9) 

Det[G h (z*)} =^(-l) p e v (^)0i 2 p(^)-©i Jvr (^J' ( 10 ) 
p 

where p is the order of a given permutations P of the set of numbers {1, 2, 3, N e }, defined by P{1, 2, 3, N e } — 
{l p , 2 P , 3 P , A P }. The single particle orbitals employed to construct the determinant have the expressions 

z*){x K {z*))\ (11) 

- ( |)a 1 + (|)a 2 , (12) 
N 

+ 1,...,0,...,--1, (13) 

A , . 

+ l,...,0,...,y-l, (14) 

and x = —irq(q — 1) /a. Note that the quantum numbers 1 (or their equivalent k) of the set of one particle functions 
01 take N e = N 2 /q 2 values. That is, the subspace formed by this set of functions is filled in the proposed may particle 



<dk(z*) = exp(- 


i h 


1 = r 2 n x 


k = 


N 


N 


?i = -yr 




N 


N 


?2 = -yr 
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fraction of the ansatz state is just v = l/q. The functions x entering the above definition are defined as follows, 



state. Since the total number of flux quanta traversing the unit cell of the periodicity region is $ = q the filling 



Xq (z*) = exp(z^i z*) x llM^i** -R*- C q )l - r*), (15) 

R 

n a i * \ , q T * a na\ 

C q = — (q.air -q.a 2 ) + , (16) 

where the vectors R (note that the vectors will be designed either by their boldface symbols of by usual symbols 
associated to their complex representations: See Appendix A) defining the product, and q are given by 

7T 

R = r\a + r 2 a 2 , R* = r\a + r 2 a 2 , a 2 — aexp(i— ), (17) 

o 

N N 

rx,r 2 = --,...,Q,...,--l, (18) 

q= A Sl + A S2 ' (19) 

N N N 
qi,q 2 = --,-— + l,...,0,...,j-l, (20) 

a 

and t?i is the elliptic Theta function 

oc 

2 ^(-1) B exp(i7ra(n - -) 2 ) sin((2n - l)u), (21) 

n=l 

Im(a) > 0. (22) 

The above functions x are representations in terms of the Elliptic Theta functions of the composite fermion like 
states derived from the exact HF single particle solutions discussed in the Ref 11 ' 12 i 13 ' 1 i 15 ' 16 , after factorizing and 
canceling zeros situated at momentum independent positions. Note that the momenta argument of the functions x m 
he definition of the orbitals in pip has been divided by q. this selection was heuristically motivated by considering the 
power q of the composite fermion should be expected to determine a momentum of the single particle wavefunction 
being q times greater than the one associated to each of the q composite fermion constituent blocks. 

It should be noted, that the case of HF states associated to a crystal periodicity showing one electron per unit cell 
(7 = 1) and v = l/q was considered here. In this sense, it can be underlined that a similar discussion to the one given 
here seems possible to be done for the case (7=5) and v — l/q. 

In the next subsection we will show that the functions Oi satisfy the following closure properties under a magnetic 
translation in an arbitrary vector of the form 1': 

T{ 1', **) 0,(2*) = exp($,(l')) 0[i-i' ]red (2*), (23) 

where [1 — \'] re d means the vector equivalent to 1 — 1' modulo an element of the lattice generated by the unit cell 
vectors ai and el 2 . The satisfaction of the periodic boundary conditions in the von Karman periodicity region of the 
functions 0i will be shown in Appendix B. 

A. Action of the magnetic translations on the functions Q\(z*) 

Let us now argue that the space generated by the set of functions ©i(z*) for all 1 values of the quantum numbers 
is closed under the action of arbitrary magnetic translations in any vector of the form 1' which are the spacial shift 
defining the Hamiltonian in This action can be expressed as follows 

T(l', **) 0,(z*) = exp(-^(y + z") - ix(z + I'*)) x ( X k{z* + l'*)) q . (24) 
But, the appearing power q of the functions x can be explicitly written in the form 
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( x ^z* + l'*)f = exp(^M'* + '-^z*) (Re^iz* -R*- \a\ - (f - r'))| - r*)J . (25) 



9 



After substituting this expression in (j2"4")) . it follows 

T( 1 , z ) 6i(z ) = exp( — — H — / - zxri )x 



exp(-ix **) (exp( l z^H^fyz* R* \a* 2 q(l* l*'))\ t*)J 



(26) 



I*' , ,(1) 



= exp(-^^- + - tad'*) x 0i_,/(2*). (27) 

As it can be noted, the last line in the above equation is the same analytic expression defining the function Oi, 
but in which the shifted quantum number value 1 — 1', may lay outside the set defining the allowed values of such 
quantum numbers. However, the difference 1 — 1 can be expressed in terms of a proper value of a quantum number 
modulo a linear combination of some basic quasi periods of the functions O. For seeing this property, let us note that 
I* — I'* can always be written as 

1* — I* = [1* — l*]red ~ h a ~ ^2 a 2J (28) 

where [I* — l'*] re d has the quantum number form 

[r-l%ed = %a+ ga$, (29) 
i i 
N N N 

ri ' r2 = -V-2^ + 1 '--'°'">- 1 ' 

and t\ , t 2 are defined as functions of I* and V* as the unique pair of integer numbers required to express [I* — l'*] re d 
in the above form. The definitions of 1 and 1' permits to write for them 

h = -Floor[ ^ J — , 

7 

t 2 = -Floor[ ^— ^]-, (30) 

where the operation Floor[n] as usual, is defined as the nearest integer number being lower or equal that the real 
number n. Note that the indices r 2 take arbitrarily large values because they define the translations entering in the 
periodic expansion of the Coulomb potential. On another hand the values r±, r 2 only define the set quantum numbers. 
Then, we can write for the main quantity entering in the arguments of the Theta functions as follows 

(I* - I'*) + |a* + R* = [1* - l*'] red + |aS + R* - ha - t 2 a*. (31) 

Therefore, the shifts in the vectors V* produce a new argument which coincides with an argument of another function 
in the considered set, plus a shift in an integer number of the unit cell vectors ai and & 2 . However, the functions 
01 under consideration satisfy exact recurrence relations under such displacements. They are direct consequences of 
the usual transformation properties of the Theta functions under shifts in their quasiperiods. The derivation of these 
relations is presented in Appendix B. In order to start making use of this property, the expression (|26p for the action 
of the translations on the functions 0\ can be written in the following form 



j(Y) /*' id) 

T(l', z*) 6,(z*) = cxp(-^^i- + -^-l'* - (32) 



exp(-i>r z*)exp( ' £ 2 * ^ z*) x (fi(«* (/', I, z*) + ha + t 2 a*)) q , 
tt(w*) = l[6 1 (^(w*-R*)\-T*), (33) 

R 

v*(l',l,z*)=z*-la*-[l*-l* , ] red . (34) 
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Then, we can make use of the relation (|B12[) in Appendix B, for the above defined functions f2 

v * 

Sl(v * + m*) = exp(iirm 2 + 2i irm2— + i nT*m 2 (m 2 + 1) ) Q(v*), (35) 

a 

in which the vector m for our situation should be chosen as m = rriiSLi + m 2 SL 2 being equal to ti&i + *2 a 2- Applying 
this formula for transforming the product V\ of Q functions appearing in the previous expression: 

Vi = (Q(v*(l', I, z*) + ha + t 2 a* 2 )) q (36) 



we can write the shifted f2 functions as multiplicative factors of the original ones as follows 

Q(v* (/', l, z*) + ha + t 2 a* 2 ) = exp(?;7T< 2 + 2wrf 2 " + wri 2 (l + *2))0(u* (/', Z, z*)). (37) 



Then, Pi have the expression 



Pi = exp(i7T^ 2 + 2wrgt 2 - 1 ' ' Z ' + inqt 2 (l + t 2 ))x (38) 



((n( v *(i',i,z*))) q . 

Therefore, the action of the translations on the function Oi can be expressed as 



T-l/l/ *\ r\ / • ^ I • CVS ^ ,/:! 

i (1 , z ) <o\{z ) = exp( — 1 — — j h « «>ft — sxz ) 



I 1 m 1/1 s * * c a l a 2 n 
exp 1- 1 \red)2 z 9^1 1" *2 J2 Z 

\ r o K q q 

exp(iirqt 2 + -^^1{ Z * - ^ p-/']red) + i^T*qt 2 {\ + t 2 ))x 

a 2 



|a|-[r-H rerf )|-T*)J . (39) 
After employing the following identity 

(a 2 ) 2 a = q X 2^, 

which reflects that exactly g flux quanta pass through the unit cell formed by the vectors ai and a2 , by also noticing 
that (i~) 2 = 0, because the vector ai = a(l,0, 0), the following closure relation arises 



(in ;/* (U 11 

T(l',z*) 61(2*) = ^{-1-4^ + 1—9 ^'*)x 

2r„ ^ z 



2iriqt 2 qar* . . 

exp{tTTqt 2 (— h [l-l \red) + 9*2(1 + t 2 ))x 

a 2 

exp(-ix**) exp(-^([l - \'} red ) 2 z*) (^[O^z* - R* - |a* - q[l* - l*'] red )\ - t*)^ 

= exp($ 1 (l'))e [1 _ r]rei ( 2; *). (40) 
The space independent multiplicative factor has the formula 

exp($ 1 (l')) = expH^L + ixl '*) x 

exp(mqt 2 - ^^(^- + [l-l'] red ) + iTTT*qt 2 (l + t 2 )). (41) 
a 2 

The integer numbers h and t 2 were fixed before in (f3"0"l) as functions of 1, 1' by the condition of enforcing the vector 
[1 — l']red to be an allowed quantum number. 
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B. The eigenstates of the FQHE in the LLL for projected Coulomb interaction 

Let us consider now the central point in this work. Note that the action of an arbitrary translation entering in the 
definition ([3]) of the periodic potential, on any of the functions Q\, reduces to a permutation of the set of quantum 
numbers times a space independent exponential factor. Thus, the defined subspace of functions on the LLL remains 
invariant under the class of translations appearing in the expansion of the interaction potential. 

Let us recall here that the number of inequivalent translations in vectors V is N 2 /q 2 , which is q times smaller 
than the number of flux quanta q(N 2 /q 2 ) passing though one von Karman periodicity cell of area n.Li xL^g 2 . By 
inequivalent translations we mean the whole set of translations obtained for arbitrary values of the quantum number 
1. Note that due to the closure relation, the translations differing from tho e inequivalent ones in an arbitrary 
linear combination with a integral coefficients of the vectors ai and a 2 , reduce to inequivalent translation times 
a multiplicative factor. For q = 1 the set of translations of the form I' transforms the LLL level in itself since the 
number of translations coincides with the number of independent function in this level and any translation implements 
a permutation of the functions. However, for q > 1, the number of translations I' is again equal to the number of 
functions Q\, but this common number is q times smaller than number of independent functions in the LLL (equal 
to the number of flux quanta passing the region of area n.Li xL 2 /<7 2 ). 

Let us now inspect how the Hamiltonian of the system acts over the wavefunction \& defined in @ . Through the 
help of relation ([23]) the result of this action can be written in the form 

i <2irp 2 r 2 1 1* 

H *(zl,z* 2 , z* N J = — £ ex P(~^) *?) T(-l, z*)x 

c i<j 1 ' ' ° 

£(-ire liP {z{)& hP (z 2 *)...e ljvf (z* N j, 
p 

1 2ne 2 r 2 , 11*, 

i<j 1 1 1 ° 

J2(-l) p 6 V (zt)Q l2P (z 2 *)... exp($ liP (l))6 liP _ 1 «)... 
p 

exp($ I . p (-l))e 1 . p+1 (z*)...e IjVe (^ e ). (42) 

Let us now recall that the action of the Hamiltonian on an antisymmetric function of the coordinates gives again 
an antisymmetric function of the coordinates. Therefore, it should be invariant under the operation of adding its 
expression for all permutations of the arguments z*, z 2 , ...z^ , after multiplication by the order of the permutations, 
and a division by iV e !. However, it can be noted that the action of the two translations involved produce other functions 
of the same set, showing modified momenta quantum numbers and being multiplied by constant exponential factors. 
Therefore, when the increment 1 is unable to transform both quantum numbers IjP and l^p between themselves modulo 
a vector of the form mi&± + rn 2 a-2 {mi and m 2 integers), the corresponding antisymmetrized term in the sum over 
permutations in equation (I42|) will vanish. This is a result of the fact that, in the mentioned situation, there will be 
two coincident states in the Slater determinant resulting from the above described auxiliary symmetrization process. 
Therefore, these properties imply the following eigenvalue equation to hold 

H ^ = E P V , (43) 

i<j IP ' ° 

exp^W (1) + <&« (-1)) x S^ p \l (Lp - L,p)), (44) 

where fr( K > p ) (1) = 1 if 1 = modulo a lattice vector of the form miai + m 2 a 2 , and vanishes otherwise. Equation 
([4"3"| is the main result of the present work. It states that $ defined above is an exact eigenfunction of the many 
electron system in FQHE regime for a periodically extended Coulomb interaction. Note that the solution is obtained 
for any finite particle systems satisfying the conditions assumed here. Thus, it can be employed to investigate the 
thermodynamic limit. The evaluation of the energy per particle associated to the obtained solution will be pursued 
in paper II. 

Let us briefly comment here about how the large number of translations entering the definition of the periodic 
potential in the growing number of particles limit, is however able produce the closure relation which was central in 
obtaining the argued result. 



8 



It can be noted that for q = 1, the set of translations of the form 1' transforms the LLL level set of functions satisfying 
the assumed periodic conditions, in itself. That is, the number of translations coincides with the number of independent 
function in this level which satisfy the boundary conditions, and any translation implements a permutation of the 
elements of a special complete basis of eigenfunctions </> q defined in Refsj 10 ! 11 ' 31 . These orbitals are eigenfunctions 
of the magnetic translations on a lattice showing exactly one unit of flux per unit cell. The form of this unit cell 
can be arbitrary, the only requirement on it, is to being pierced by one flux quanta. An important property of those 
translations which is argued by example in Refsj n i 31 , is that any magnetic translation in the LLL is equivalent to a 
shift in the momentum quantum numbers of the functions in the mentioned basis. Let us consider such a one flux 
quantum unit cell as formed by the vectors ai jq and a 2 . The number of such cells in the area of the sample we are 
considering is exactly q(N 2 /q 2 ), that is, the number of flux quanta passing through the sample area. In this case, 
the momenta associated to the complete set of functions (/> q in the chosen sample area ar o 10 ' 11 ! 31 



q = 


9i . 92 

9 






Sl = 


n x a 2 

2 ' s 2 = 

7 o 


n 


x a 2 

^ o 


qi = 


N TV 

-T'-T + 1 - 


,0,. 


TV 

'T 


92 = 


TV TV 
~Yq^Yq + '■ 


..,0, 


TV 



(45) 

Let us consider now the (-y) 2 translations generated by the vectors 1' when 1' is restricted to be a quantum number 
of the functions Then, the shifts in the momenta q produced on the functions q after acting on them with a 
translation 1' is given byii 

5q=-^A(l') 

he 

= n x (1 ) 

he 2 V ' 

= 4n x 1' 



f o 

rj TV N 
= (-QjjZi + Js 2 ). (46) 

9 

This expression indicates that when q > 1 there exist a proper subset of the functions </> q that are closed under all 
the magnetic translations defining the specially periodic Coulomb operator. This conclusion follows form the fact 
that the shifts (5q are closed under addition, but they only have (^) 2 values within the first Brillouin zone, of the 

total number q{^) 2 of momenta forming this zone. Therefore, the existence of a close proper subset formed by the 
functions Q\ within the LLL subject to the considered periodic conditions is allowed. 



C. Formula for the energy per particle of the exact eigenstates 



Let us employ in what follows the periodicity conditions for finding a formula for the energy per particle in terms 
of the 9i functions, for the FQHE eigenstates. Since the state ^ is a Slater determinant the total energy in the von 
Karman periodicity zone can be expressed in the form 

E = < * | V c \ * > 

= f f d/i(x)d/ i (x')V r c (x-x / )0(x,x / ), (47) 
Ja c j a c 

5«e( X > X ') = P( X , X )P( X ', X ') - P( X , X 'M X '> X )) 

P (x,x') = ^e 1 (x)er(x'), (48) 

1 
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where 0i(x) as a function of the coordinate vector x is naturally defined by 

9i(x) = 0i(((x)i, (x) 2 )) = 6,(z*) = ei((x)! - i(x) 2 ), (49) 
and as usual, in the Landau gauge the integration measure is given by 

d/*(x) = d x exp(- ^ x y ). (50) 

We also assumed here that the functions 0i(z*) were already normalized in the von Karman periodicity zone. 
Let us consider in what follows the following type of integrals 

1= f d/i(x)K c (x-x / )0i(x)e 1 /(x), (51) 
J A c 



which are building elements of the expression for the energy (|47p . To start, note that the region of integration is 
limited to a cell of the lattice R 9 , which was named above as A q . For this purpose, we can rewrite (|B5[) . after 
substituting the explicit form of the magnetic translation invariance in the Landau gauge (|A6j) . as follows 

Q 1 (z*) = T(K q ,z*)e 1 (z*) 

,(R q ) 2 R q .(R q ) 2 z* . 
= exp(-z — ±-i — H — — )6i(z +R 9 ). (52) 

Employing this relation for the two functions 0i and Q\i in ([51) and taken into account the periodicity of the 
potential in the lattice R q , it follows 

I=( dM(x)F c (x-x')ei(x + R g )6r(x + R 9 )x 

JA C 

jR q ) 2 (R* q - R q ) , (R q ) 2 (z*-z) , 
exp(-« r-| % — 2 

= [ d M (x + R ? )F c (x-x')ei(x + Rg)6v(x + R 9 ). (53) 

This expression explicitly shows that the integral I performed over an arbitrary cell of the lattice R q , give the 
same result when performed over any of such cells. This property directly check that total energy can be evaluated 
as an integer limited to the first cell of the lattice R q as the von Karman conditions implies. Then, the total energy 
has the formula 

E„_ = ( ¥ I V c \ * ) 

d/i(x)cfyi(x')V c (x - x / )sf JVe (x,x / ). (54) 



a c JA,_ 

The evaluation of the energy per particles of the finite systems, their thermodynamic limits and the pair correlation 
functions will be considered in paper II, which is now in preparation. 

D. Remarks on the thermodynamic limit 

Let us briefly comment in this subsection on the thermodynamic limit of the determined exact wavefunctions. 
As known, in the large magnetic field regime in which the system is projected in LLL level, the only term in the 
Hamiltonian is the Coulomb interaction. Here we assumed von Karman periodic conditions in a magnetic field 
(See (H), with spacial periods Li/g and L 2 /q . Let us assume as usual , that the systems has a homogeneously 
distributed jellium which compensates the net charge of all the particles contained in it. In this situation, when 
the |Li|/<7 = | L2 1 / <? = L/q = aN/q tends to be large in the thermodynamic limit, the interaction between any 
pair of particles sitting at arbitrary but finite distances from the origin, should approximate the usual Coulomb 
potential, as -| becomes very much greater than the distance between the mentioned two particles. This property is 
expected to be valid because when L/q is very much larger than the distance between the particles, their interaction 
potential becomes negligible assumed the compensating fields of the jellium homogeneous charges are considered (or 
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equivalently, imposing that the inter particle potential is assumed as short ranged, by substituting the Coulomb one 
by a Yukawa screened potential at large distances). In other words, the jellium charge, should reduce the fields of 
the far away laying image charges to be dipolar or higher multipole short ranged contributions, which should produce 
vanishing forces in the thermodynamic limit. Therefore, it can be expected the solution found in last section, when 
considered in the limit of a large number of particles, should describe the FQHE systems when the limit L/q — > oo 
is taken. In this case of large samples the following formula for the energy per particle arises from expression (|54| for 
the total energy of arbitrary finite systems 



Exact eigenstates of a two dimensional electron gas (2DEG) in a magnetic field are constructed for finite samples 
of arbitrary sizes, and periodic extensions of the Coulomb potential after restricted to the LLL level. Magnetic von 
Karman periodic boundary conditions are imposed in the Landau gauge. Formulae for the energy per particle and the 
two particle density matrix of the realistic FQHE state in the thermodynamic limit are also determined by considering 
the large limit of finite systems. The found states are associated to filling factors v = l/q, for odd q. The results 
suggests the presence of integrability properties in FQHE problems. The energy per particle and pair correlation 
functions of the states for the finite system and the thermodynamic limit ones, will be evaluated in a second work 
(paper II) which is now in preparation. In the case that the energies per particle associated to the wavefunctions 
result to be lower than the estimated values for the Laughlin variational states, they could surprisingly furnish the 
ground states of the FQHE problem. However, in the opposite outcome, the determined wavefunctions could be also 
helpful in describing experimentally detected phases shown by FQHE samples at variable temperature experiences. In 
any situation, the obtained kind of solution constitutes a valuable primer of exactly solvable model in the context of 
the FQHE theory. The solutions could also be exact realizations of the early proposed cooperative rings of exchange 
states^. The evaluations of the energy per particle and pair correlation functions of the identified eigenstates will 
be presented in a related work (paper II) which preparation is in progress. 

The special one particle orbitals constructed here, were also employed (See Appendix C) to construct ansatz wave 
functions for a class of low energy states having a similar structure as the Laughlin's ones. That is, showing zeros of 
the form (zi — Zj) q when the coordinates of the particles coincide. Then, these states exhibit a similar sort of short 
range behavior as the Laughlin states. However, since the basic elements in their definition are the HF one particle 
orbitals, we expect that long range correlations incorporated in the original HF solution can remain in the proposed 
ansatz states, helping in this way to slightly lower the energy per particle. 
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APPENDIX A: BASIC DEFINITIONS 



1. 



General definitions 



In this work, bold symbols represent vectors with components v =((v) 1 , (v) 2 , (v) 3 ) and 2-dimcnsional vectors are 
equivalently defined by the two complex numbers 



v = (v) x + i(v) 2 , v* = (v) 1 -i(v) 2 . 
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Let us consider the plane containing the 2DEG. One kind of useful periodicity planar box employed in previous works 
consists in a parallelogram of sides Li = N a.\ and L 2 = N a 2 , where the unit cell vectors a! and a 2 arc defined by 

ai =0(1,0,0), a 2 =a(i,^,0), (Al) 



n= (0,0,1), a = y^r c , (A2) 

where q is an odd integer, n is the unit vector normal to the 2DEG plane, and r =^j^ is the magnetic length. 

The Landau gauge A = B{— a; 2 , 0, 0) is employed for the vector potential in accordance with Ref. [3] • Note that the 
cell with sides ai and a 2 intercepts exactly q flux quanta of the magnetic field B = _B(0,0, 1). Also, N 2 particles 
laying inside the above defined particular planar box correspond to a filling factor v = l/q. Note that this box will 
not be the von Karman periodicity box considered for determining the exact solutions, which will have sizes given by 
lix/q and L 2 /<7. The electric charge e is assumed to be negative. For convenience, N will be assumed to be an even 
positive integer number and also an exact multiple of q. Is helpful to define a spatial lattice given by the vectors 

R= i?iai +-R 2 ai, (A3) 

in which i?i and i? 2 are arbitrary positive or negative integers. The reciprocal lattice cell vectors associated to this 
lattice has unit vectors 

n x a 2 n x ai 

si = 5-, s 2 = 5-, (A4) 

qr z Q qr z Q 

Sj.aj = 2TrSij, i,j = 1,2. (A5) 
The free Hamiltonian of the single planar electron in the Landau gauge will be 

where the generators of the magnetic translation symmetry operations are 

„ ie . ieB 
G = V-— A+-— n x x, 

he he 



[G,H] = 0, 

ieB 
He 



( G )i.( G )j] = l-e y *( n )» 



The structure of the free electron wavefunctions of the LLL in the Landau gauge is 

y 2 

<p(x) = <f>(z*) exp(-^), y = (x) 2 

2 =( x )l + ( x ) 2 ^ Z * = ( X )l - ( X )2«- 

where 4>(z*) is an analytic function of z*. In the axial gauge, the electron states and the vector potential are 

x 2 

<p A (x) = <p A (z*) exp(- — ), 

l ±T 

A A = - B n x x, 

with (/)a(z*) an analytic function of z* different from <p^(z*). The general connection between wavefunctions in both 
gauges, and between the analytic forms associated with the LLL are 

ieB 

ip A {-x) = exp( — (x) 1 (x) 2 )<^(x), 
A (O=«p(^)^*). 
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Then, the Landau gauge versions of the Hamiltonian H, and the symmetry translations generators G, can be obtained 
from the similarity transformations 

h a = exp(^-(x) 1 (x) 2 ) J H"exp(-^-(x) 1 (x) 2 ) 

= — p--a a 2 = -— 

2m c 2m 

GA = CXP (5k : ( X )lW2) GcX P(-5-( X )l( X )2) 

_ ie 
= V+— A^, 

_ „ ie . 

n A = v-— a a . 

he 

These quantities have the following properties 



[(G A ) l ,(n A ) J } = o 

[(G A )i,(G A )j] = -^e U j{Ta) u 

' O 

[(U A )i, (n A )j] = -^e ai (n)i 
[H, G A ] = 

T A (L,x)=exp(— (x) 1 (x) 2 )T(L,x)cxp(- — (x)^^,). 

Helpful quantities are also the lowering and rising operators U± and G± defined by and satisfying the following 
relations 

n± = {u A \ ± (n A ) 2 i, 

[ir,n±] = ±-^ T n±, [fr,G±] = o, 
[n +) n_] = l, [G +> G_] = -i 

They have the explicit forms in terms of the complex representation of the coordinates 

Zj/ o o 

Now let us consider a property that becomes useful for projecting the Coulomb interaction in the LLL, since it is 
convenient for studying the high magnetic field limit. Firstly, let us express the scalar product of a 2D coordinate 
vector x with a vector Q, as 

tc 

iQ.x = -21.— Aa(x), 
he 

1 = r 2 n x Q. 

But the gradient and the axial gauge vector potential can be expressed in terms of the magnetic translation (MT) 
symmetry and the operator II ^ as 

\7 = l(U A + G A ), 

IP 1 

-A,(x) = -(G A -n A ). 
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Then, the exponential exp(iQ.x) can be represented as follows, 



exp(iQ.x) = exp(-21.— Aa(x)) 



exp(l.n j4 ) cxp(-l.G A ) 

ex P (i(rn+ + ra_) cx P (-i(PG+ + zg_) 



II* I 1 

CX P(~^2 ) ex P(2 rn +) ex P^ 



exp(-i(Z*G+ + /G_)), 



where the well known identity 



exp(A + B) = exp(A)exp(B)exp(-^[A,S]), 

valid when [A, B] commutes with A and with B, has been used. Applying the one particle projection operator in the 
LLL: P = J2 n \ n )( n \ ({l n )} being any complete set in the LLL), it follows that 

P« exp(;Q. Xi )pW - exp(-^)exp(-l.G A ) 

= exp(-^)exp(-l(rG + +/G_)) 
= exp(-^)exp(-i(n2fl,;+^)+ 

2. LLL projection of the Coulomb interaction operator 

Let us next consider the Fourier expansion of the 2D function /(x) which is periodic in a parallelepiped with sides 
Lf and Lf 



/(Q) = y dx /(x)exp(-iQ.x), 
5 P (x - x') = J- cxp(*Q.(x - x')), 



iceM g 

A ce;/ = n.Lf x L|. 



It should noted that here in this subsection A ce u, Lf , and are the area and the unit cell vectors of a general periodic 
lattice. Employing this expansion, the Coulomb interaction Hamiltonian can be expressed as follows, 



v °= E Eu- 



e 2 



/, . l Xj X J + Rc ' 

i,j=l,2,...N e 

1 \ ^ \ ^ 27re 2 , 

XT E E ToT ex P( J Q-( x i- x j))- 
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where R c is the lattice generated by the unit cell vectors Lf and Lf and Q are Fourier components of the periodic 
potential in the box with area A ce u. But applying on Vc the many-particle projection operator in the LLL 



p (xi,x 2 ,...,x JV j = n p o ( 



N e 

) 



i=l 

the projected Coulomb interaction Vc,a can be written in the axial gauge as 

V c ,a = PoVcPo = ^^ (x i )exp(iQ.x i )P (x i )P (x j )exp(-iQ.x;.)P (x J ), 

CeU i<3 Q 
ce(Z 4<j . ,11 o 

1 

ex P (-i(r(2^ + ^) + Z(2^.-^))), 
l=r^n x Q. 

This operator may be written in the representation given by the analytic factors by considering the following rules 
expressing how the translation operators transform when extracting the exponential factors from the wavefunctions: 

Then, the projected Coulomb interaction in the axial gauge acquires the form 

K.J 1 

ex P (-ra. + (i) 



i<3 1 

W )=exp(-^-g)e*p(X.A), 

where Ta(1, z*) is the translation operator in the axial gauge acting in the space of analytic functions. After 
performing the many-particle similarity transformation which transforms axial gauge operators into the Landau 
gauge, the Coulomb interaction reduced to the LLL in this latter gauge takes the form 



yLan.au = ( cxp( _ L^_ } } yAnal ( JJ ^(Li^)) 



AEEt exp(-^)r(-ix)r(i,,;; 

Acell i<} j |1| ^ Q 



where now the Landau gauge translation operators acting on the analytic parts of the single particle wavefunctions 
have the form 

T(l, *•) = e*p(-i^f - iQhg.) exp(r A). (A6) 
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APPENDIX B: VON KARMAN BOUNDARY CONDITIONS 

Let us argue below that the functions 0i(z*) defined above are invariant under the translations. 

T(^,z*)=T(N^,z*). 

q q 

This property, implies that those functions satisfy von Karman periodicity condition under translations in the periods 
^-.This result can be directly derived from the relation 

Consider first the shift in Then, let us fix the parameters in (|40|) in the form 

r-±. 

q 

i-i' = i-^i. 

q 

[1 - l']red = 1, 

N 

h = -,t 2 = o, (Bl) 
q 

where it has been employed that, since Li/g=JVai/q is a linear combination with integer coefficients of the vectors 
a! and a 2 thanks to the fact that N was supposed as even and also as a multiple of q. Substituting these values in 
(|40| directly leads to the first von Karman boundary condition 

T(^,z*) 0!(z*)= Qi(z*). (B2) 

Similarly for a shift in ^p, the required parameters in (|40|) become 

q 

1-1=1- i^, 

q 

[1 — 1 ]red = I) 

N 

h=0,t 2 = -, (B3) 

q 

which after performing a little more complicate algebra implies the second von Karman boundary condition 

T(^,z*) ei(z*)= Qi(z*). (B4) 

q 

Therefore, the functions 0i are invariant under magnetic translations in the lattice vectors pertaining to R q as 

T(R q ,z*)e 1 (z*) = e l (z*). (B5) 

1. An auxiliary relation 

Let us consider the product of Theta functions 



R 

where to* is the complex representation of the vector m = mi a!+TO 2 a 2 , and the product over R runs as before over 
the vectors 
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R = i?ia!+i? 2 a 2 , (B6) 
N N 

Rl,R 2 =-y, .-.,-1,0,1, ...y-l- (B7) 

The simple properties of the Theta functions under shifts in vectors of the form RiSlx, and the even character of N, 
make independent of R\. Therefore 



+ m*) = Q(v* + TO2I2) 

= - {Ria + (R 2 - m 2 )a*))| - r*). (B8) 

R 

Now, it can be noted that 7712 can be interpreted as a shift in the arguments R of the factors, defining a modified set 
of arguments R'. Then, let us write the product V of all the factors in which the modified arguments R! are not 
appearing in the original set R. Assuming 7712 > 0, this product can be given the form 



"12 nr 

fli 3=1 

= n n - ( R i a + (y - jk)) + ttt-i - r*). 

fli 3=1 

After iterating the usual formula of the 0\ Theta function under a single shift in 7tt* , the transformation rule when 
the argument is shifted an integer k times 7rr*, can be written in the form2S 



(~l) fc 
q k2 (r*) 

q(r*) = cxp( — Z7TT* ) , k = —00, —1, 0, 1, ...00. 



n(u* + knT*\ -r*) = \ 2 ' exp(2iku*)6 1 {u*\ -t*))0i(u*1 -r*) (BIO) 



The above expression allows to write V in the form 



in 2 



v = n n t4 ^t^* - (i?ia + ( f ~ x 

Ri 3=1 q{T ' 



-(fl l0 +(^-jK))|-r*) 



L v v 2 

f* 2iirTn 2 
cxp(+2i7rm 2 - — — > itio- 



ma ma ^ 

2z7tt* e ^ x n n - (i?ia + ( y ~ ~ r * } - 

3=1 -Ri 3=1 



Then, the use of the relations 



— — 1 
2 x 

Rl=-f 

3=1 



"T- 2 (m 2 + 1) 



reduces the expression for to the same product in which the mi shift is absent, after multiplication by an exponential 
factor, as follows 
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V = II II 6 ^ v * ~ ( R i a - - jam - r*) (Bll) 

fli 3=1 

= exp(i7TTO 2 + 2j 7rTO 2 — + i 7T T* JT12 (TO.2 + 1 ) ) x 

n n °^ v * + (y - - o- 

fli 3=1 

Then, the following resulting relation arises 

f2(v* + to*) = exp(i7TTO 2 + 2i Trm 2 — + i 7tt*TO2(to,2 + 1) ) fl(v *), (B12) 

a 

checking that shifts in the vectors to* of the functions f2 reproduce the same function times an exponential factor. 

APPENDIX C: ANSATZ STATES FOR THE FQHE PROBLEM SHOWING TRANSLATION 

SYMMETRY BREAKING AT v — - 

i 

Let us employ in this section the composite fermion-like basis defined in Refs . 14 ' 15 ^ 16 and in (fTBf . to construct a 
special class of many-particle state associated with the fractional filling factor l/q. These wavefunctions show a two 
particle pair correlation function having the known optimized behavior g(xi — X2) ~ c(z* — z^) q which is helpful in 
reducing the correlation energy. However, the original source of the definitions for the basis functions Xk ; (which 
come from HF single particle solutions after extracting the fixed position zeros) leads to expecting the presence of 
long range correlations further diminishing the energy per particl o 14 i 15 ' 16 . The ansatz states will be simply defined 
by 

y G {.z{,z* 2 ,...z Nc )=Det[Q h {z* j )}< 1 , 

Det[e u (z*)} = ^(-i)fe liP (z*)e l2P (z 2 *)...e ljVeP (z^) 

p 

These wavefunctions satisfy the boundary conditions 

T(L 1 / g ,z*)* G (z 1 *,z*,...z^J = ^ G {zl,zl,...z* Nc ), 
T(L 2 /q,z*)'f G (z* 1 ,z*,...z* N J = * G (z*,z*,...z^J, i=l,...,JV ei 

which directly follow from the boundary conditions for the functions Oj. 



1 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 



2 



B. I. Halperin, Z. Tesanovic, and F. Axel, Phys. Rev. Lett. 57, 922 (1986). 
3 S. M. Girvin and A. H. MacDonald, Phys. Rev. Lett. 58, 1252 (1987). 
F. D. M. Haldane and E. H. Rezayi, Phys. Rev. B31, 2529 (1985). 
N. Shibata and D. Yoshioka, Jour. Phys. Soc. Japan 72, 664 (2003). 
S. Kivelson, C. Kallin, D. P. Arovas, and J. R. Schrieffer, Phys. Rev. B36, 1620 (1987). 
J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 
F. Claro, Solid State Commun. 53, 27 (1985). 
F. Claro, Phys. Rev. B35, 7980 (1987). 
A. Cabo, Phys. Lett. A171, 90 (1992) 
A. Cabo and F. Claro, Phys. Rev. B70, 235320 (2004). 
A. Cabo, F. Claro, A. Perez and J. Maze, Phys. Rev. B76, 075308 (2007) 
A. Cabo, F. Claro and A. Perez, Phys .Rev. B66, 035326 (2002). 
A. Cabo, F. Claro, arXiv: | cond-mat/0702251 (2007). 

15 A. Cabo, F. Claro, arXiv: cond-mat/07 (2007). 

16 A. Cabo, F. Claro and D. Martmez-Pedrera, ICTP Preprint: IC/2007/092(2007). Available at: |http: / /publications.ictplt 



4 



18 



17 Y. P. Chen, R. M. Lewis, L. W. Engel, D. C. Tsui, P. D. Ye, Z. H. Wang, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 
93, 206805 (2004). 

18 G. A. Csathy, D. C. Tsui, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 92, 256804 (2004). 

19 W. Pan, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 88, 176802 (2002). 

20 Yuan P. Li, T. Sajoto, L. W. Engel, D. C. Tsui, and M. Shayegan, Phys. Rev. B47, 9933 (1993). 

21 E. J. Bergholtz and A. Karlhede, J. Stat. Mech.: Theory Exp., L04001, (2006). 

22 T. H. Hansson, C. C. Chang, J. K. Jain, and S. Viefers, Phys. Rev. Lett. 98, 076801 (2007). 

23 T. H. Hansson, C.-C. Chang, J. K. Jain, and S. Viefers, Phys. Rev. B76, 075347 (2007). 

24 E. J. Bergholtz, T. H. Hansson, M. Hermanns, and A. Karlhede, Phys. Rev. Lett. 99, 256803 (2007). 

25 E. J. Bergholtz, T. H. Hansson, M. Hermanns, A. Karlhede, and S. F. Viefers. larXiv:0712. 38481 (unpublished). 

26 G. Fano, F. Ortolani and E. Colombo, Phys. Rev. B34, 2670 (1986). 

27 S.T. Chui, T. M. Hakim and K. B. Ma, Phys. Rev. B33, 7110 (1986) 

28 P. K. Lam and S. M. Girvin, Phys. Rev. B30, 473 (1984). 

29 D. Levesque, J. J Weis and A. H.MacDonald, Phys. Rev.B 30, 1056 (1984). 

30 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (1964) Dover Publications, New York. 

31 R. Ferrari, Phys. Rev. B42, 4598 (1990). 

32 E. Anisimovas, M. B. Tavernier and F. M. Peeters, Phys. Rev. B77, 045327 (2008). 

33 M. B. Tavernier, E. Anisimovas and F. M. Peeters, Phys. Rev.B B74 125305 (2006). 

34 T. Stopa, B. Szafran, M. B. Tavernier and F. M. Peeters, Phys. Rev.B B73, 075315 (2006). 



